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Abstract 

We give a short description of the present situation of lattice QCD simulations. We then focus on the computation of 
the anomalous magnetic moment of the muon using lattice techniques. We demonstrate that by employing improved 
observables for the muon anomalous magnetic moment, a significant reduction of the lattice error can be obtained. 
This provides a promising scenario that the accuracy of lattice calculations can match the experimental errors. 
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1. introduction 

The interaction between quarks at large distances be- 
comes strong such that analytical methods as perturba- 
tion theory fail to analyze QCD. A method to never- 
theless tackle the problem is to formulate QCD on a 
4-dimensional, euclidean space-time grid. This setup 
first of all allows for a rigorous definition of QCD and 
leads to fundamental theoretical and conceptual investi- 
gations. On the other hand, the lattice approach enables 
theorists to perform large scale numerical simulations. 

In the past, lattice physicists had to work with a num- 
ber of limitations when performing numerical simula- 
tions which turn out to be extremely expensive, leading 
to the need for Petaflop computing and even beyond, a 
regime of computing power we just reach today. There- 
fore, for a long time the sea quarks were treated as in- 
finitely heavy, indeed a crude approximation given that 
the up and down quarks have masses of only O(MeV). 
In a next step, only the lightest quark doublet, the up and 
down quarks, were taken into consideration, although 
their mass values as used in the simulation had been un- 
physically large. 



Nowadays, besides the up and down quarks, also the 
strange quark is included in the simulations. In addi- 
tion, these simulations are performed in almost physi- 
cal conditions, having the quark masses close to their 
physical values, large lattices with about 3fm linear ex- 
tent and small values of the lattice spacing such that a 
controlled continuum limit can be performed. The situ- 
ation of present days simulation landscape is illustrated 
in fig. [T] taken from Q). In the figure, the black cross 
indicates the physical point. 

The fact that presently simulations close to the phys- 
ical situation can be performed is due to three main de- 
velopments: z) algorithmic breakthroughs which gave a 
substantial factor 0(> 10) of improvement; ii) machine 
development with a computing power of the present 
BG/P systems which is even outperforming Moore's 
law, Hi) conceptual developments, such as the use of 
improved actions which reduce lattice artefacts and the 
development of non-perturbative renormalization. 

As an example of physical results we can achieve 
presently, we show in fig. [2] the continuum extrapolated 
strange baryon spectrum as obtained by the European 
Twisted Mass Collaboration (ETMC) (3) of which the 
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Figure 1 : The values of the lattice spacing a and pion masses m n as 
employed presently in typical QCD simulations by various collabo- 
rations as listed in the legend. The cross denotes the physical point. 
The figure is taken from |T| where also further details and references 
to the various collaborations can be found 
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Figure 2: The continuum strange baryon spectrum from the ETM 
collaboration | 2 | using Nf = 2 flavours of quarks with only mass- 
degenerate up and down quark masses. 



Clearly, this is a very interesting result. It means that 
either in the theoretical calculation something has been 
neglected or has not properly been included. Or, some- 
what much more exciting, the discrepancy points to a 
breakdown of the standard model of particle interactions 
and the inconsistency stems from effects of some yet un- 
known new physics beyond the standard model. 

Indeed, calculations show that these new physics ef- 
fects would lead to a correction to the anomalous mag- 
netic moment of size 

Here Wi epton is the mass of one of the leptons and 
A^newphysics represents the mass (or scale) of a parti- 
cle originating from the (unknown) new physics be- 
yond the standard model. The formula in eq. ([T]i shows 
that in the case of the muon anomalous magnetic mo- 
ment the effect of new physics would show up about 
(m^/me) 2 ~ 4 ■ 10 4 times stronger than in the case of the 
electron. In principle, the r-lepton would be even more 
suitable to detect these new physics effects, but unfortu- 
nately due to the very short lifetime of the r lepton the 
experimental measurements of the anomalous magnetic 
moment of the r are presently much too imprecise to un- 
veil a possible new physics contribution. This leaves us 
then with the muon anomalous magnetic moment as the 
ideal place to look for new physics and indeed a large 
number of works has been devoted to explore this pos- 
sibility, see 1 5 1. 



authors are members. The first complete calculation of 
the baryon spectrum was achieved by the BMW col- 
laboration 1 3 1 and nowadays a number of lattice groups 
are providing calculations of the hadron masses and 
nucleon structure, see e.g. [4| for a recent review. 
The baryon spectrum calculation has been considered 
a benchmark study for lattice QCD for a long time. It 
is therefore very reassuring that finally this important 
result can be obtained precisely from ab-initio and non- 
perturbative lattice simulations. 

2. The anomalous magnetic moment of the muon 

The progress of lattice calculations discussed in the 
previous chapter motivates to address more demanding 
quantities than the baryon masses. A prime example 
is the anomalous magnetic moment of the muon = 
{gf! - 2)/2. The reason for looking at this quantity is 
that here theory and experiment disagree and one finds 
a™ 13 - fl * = 2.90(91) x 10" 9 which leads to a larger than 
3cr level discrepancy. 



3. When the lattice enters the game 

It has been found that electromagnetic and weak in- 
teraction effects can by far not serve as an explanation 
of this discrepancy [5]. However, the strong interac- 
tion can have a large effect since the hadronic contri- 
butions ajf d dominate the uncertainty of the standard 
model value. The problem is that the strong interac- 
tions of quarks and gluons in QCD are intrinsically of 
non-perturbative nature. Taking these contributions into 
account by perturbation theory is therefore rather doubt- 
ful. Employing additional model assumptions to esti- 
mate ajj ad will not provide a fully controlled and reli- 
able calculation of the hadronic contributions and hence 
an unambiguous and stringent test whether the stan- 
dard model is correct or must be extended by some new 
physics cannot be performed. 

It is exactly at this point where lattice field theory 
methods applied to quantum chromodynamics can help 
- at least in principle. In lattice QCD the theory is for- 
mulated on a discrete 4-dimensional euclidean space- 
time lattice and the theory is approached by means of 
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Figure 3: We show an earlier result for the hadronic contribution to 
the muon anomalous magnetic moment computed by the authors. 



numerical simulations. It goes beyond the scope of this 
article to explain the mathematical concepts of lattice 
QCD but such numerical simulations allow then to com- 
pute physical quantities in a fully non-perturbative fash- 
ion without relying on any model assumptions or ap- 
proximations. 

Of course, the discretization itself induces a system- 
atic error which must be removed by making the lattices 
finer and finer until the continuum of space time points 
is recovered by some suitable extrapolation process, a 
procedure which is called the continuum limit. In addi- 
tion, the simulations necessarily demand a finite num- 
ber of lattice points which can lead to finite size effects 
when the lattice is not large enough in physical units. 
Finally, often the simulations need to be performed at 
values of hadron masses that are larger than the ones 
observed in nature. The reason is that for smaller and 
smaller hadron masses the computational costs increase 
rapidly such that one is restricted to values of, say, pion 
masses that are a factor of about two larger than the ones 
observed in nature. 

All these systematic effects that appear in lattice sim- 
ulations need to be controlled in a quantitative way. For 
example, to reach physical values of the pion masses, 
an extrapolation to the physical point where the pion 
mass assumes its physical value needs to be performed. 
This appeared to be very problematic in the past. This 
is illustrated in fig. [3] The figure shows that the lattice 
results for aj) ad are significantly below the experimental 
number. An extrapolation to the physical point, recon- 
ciling the lattice data with experiment becomes in this 
situation very difficult and even needs some additional 
model assumptions. This all leads to an error of ajj ad as 
obtained from lattice simulations that is about a factor 
of 10 larger than the phenomenological one O. The 
lattice community have been therefore rather sceptical 
in the past that lattice QCD can provide a significant 
contribution to our understanding of the discrepancy in 



One additional suspicion has been that so-called dis- 
connected (singlet) contributions could be substantial. 
In all existing lattice calculations these contributions 
were neglected, however. The reason is simply that 
these contributions are very noisy and therefore hard to 
compute reliably. Nevertheless, in ref. Q a dedicated 
effort has been undertaken to calculate for the first time 
these contributions. As a result, it could be established 
that the dis-connected contributions are in fact small and 
can be safely neglected. In addition, also the effects of 
non-zero values of the lattice spacing and the finite vol- 
ume turned out to be small. Thus the difficulty to recon- 
cile lattice data with the experimental result, shown in 
fig. [5] is rather puzzling. 

A resolution of this puzzle was only given this year 
by us 0. We observed that by a suitable redefinition of 
the lattice observable needed to compute aj} ad a much 
smoother and much better controlled approach to the 
physical point can be achieved. 

To illustrate the idea, let us give the definition of ajj ad , 



£ 

Jo 



dQ 2 -^aj(r)U R (Q 2 ) . 



(2) 



Here a is the electromagnetic coupling and Ur(Q 2 ) the 
renormalized vacuum polarization function, Tl R {Q 2 ) = 
n((2 2 ) - n(0). The functional form of a>(r) is analyti- 
cally known and the argument r is given by r — Q 2 /m 2 
where m^ denotes the mass of the muon and Q a generic 
momentum. The key observation is now that on the lat- 
tice there is a large freedom to choose a definition of 
r. The only requirement is that in the limit of reach- 
ing a physical pion mass the continuum definition of 
r — Q 2 /m 2 is recovered. Hence, one may define 



flat. = Q 2 



//phys 
H 



(3) 



with possible choices for H 
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ri : 
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(4) 



Here, my(mps) is the mass of the p-meson and /V(mps) 
the p-meson decay constant as determined on the lattice 
at unphysical pion masses mps. Furthermore, m p and f p 
denote the corresponding p-meson mass and decay con- 
stant at the physical point. All the definitions in eqs. Q 
guarantee that indeed the desired definition of r is re- 
covered in the limit of a physical pion mass since then 
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Figure 4: The present day status of con^n,,.- „„ 
using the improved observables discussed in the text. The curves cor- 
respond to the definitions r\, ri and r^ in eqs. |3J from bottom to top. 



by definition m v (m P s) and /y(/«ps) assume their physi- 
cal values. In fig. H we show the results for a^ A for all 
three definitions of r. Clearly, for the definitions r2 and 
r? the behaviour of the lattice data towards physical pion 
masses is simply linear and allows for a controlled ex- 
trapolation to the physical point. As a result, one finds 
using definition r2 in eq. (|4]) values from the lattice com- 
putations and experiment 



a 



had.latt 

had,exp 
H,Nf=2 



5.72(16)- 10 



-8 



5.66 (5) • 10 



-8 



(5) 



In the equations above, the index Nf = 2 indicates that 
in the lattice QCD calculations only a mass-degenerate 
pair of up and down quarks were used. Since the strange 
and charm quark flavours were neglected the simula- 
tions do not correspond to a fully physical situation 
which also leads to some ambiguity in the experimen- 
tal extraction of a had ,' ex ~. This shortcoming needs to be 
overcome in the future. 



4. Conclusion 



In conclusion, using the modified and improved def- 
initions of a{; ad on the lattice it is not only possible to 
recover the experimental result. As the comparison in 
eq. |5]l shows it is now also possible to come signifi- 
cantly closer to the experimental accuracy. The idea of 
the improved observables which led to a much reduced 
error for an important quantity such as g^ - 2 has there- 
fore led to the promising situation that lattice QCD cal- 
culations can match the experimentally obtained errors. 

As already mentioned above, the results for ajj ad here 
discussed have been achieved for the case of two mass- 
degenerate quark flavours. What is needed in the future 
is the inclusion of the strange and the charm quarks to 



allow for a direct comparison to the experimental re- 
sults. In addition, newly planned experiments at Fer- 
milab [8] and J-PARC [9| are aiming at an accuracy of 
below 0.5% for the hadronic contribution to the muon 
anomalous magnetic moment. To match this accuracy 
dedicated lattice simulations have to be performed on 
large enough volumes and as close as possible to the 
physical point. In addition, explicit effects of isospin 
breaking and electromagnetism might need to be in- 
cluded. All this is in principle reachable within lattice 
QCD but, it constitutes a real challenge for the lattice 
community. An even larger challenge are contributions 
to gp - 2 that appear at higher order of the electromag- 
netic coupling, most notably the so-called light-by-light 
contributions. 

However, a number of lattice groups J7] [10] fTTTl are 
working on this problem presently and it can be ex- 
pected that the lattice will provide a significant contri- 
bution to answer the question, whether the observed dis- 
crepancy in the muon anomalous magnetic moment is 
indeed a sign of new physics. 
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